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Abstract 

In this paper we derive analytically the evolution equation of the interface for a model of sur- 
face growth with relaxation to the minimum (SRM) in complex networks. We were inspired by 
the disagreement between the scaling results of the steady state of the fluctuations between the 
discrete SRM model and the Edward- Wilkinson process found in scale-free networks with degree 
distribution P(k) ~ k~ x for A < 3 [Pastore y Piontti et ai, Phys. Rev. E 76, 046117 (2007)]. 
Even though for Euclidean lattices the evolution equation is linear, we find that in complex het- 
erogeneous networks non-linear terms appear due to the heterogeneity and the lack of symmetry 
of the network; they produce a logarithmic divergency of the saturation roughness with the system 
size as found by Pastore y Piontti et al. for A < 3. 
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During the last few years the study of complex networks has moved its focus from the 
study of their topology to the dynamic processes occurring on the underlying network. This 
is because many physical and dynamic processes use complex networks as substrates. Re- 
cently, many studies of dynamic processes on networks, such as epidemic spreading traffic 
flow [3, |4J , cascading failure [5( , and synchronization , have demonstrated the impor- 
tance of the topology of the substrate network in the dynamic process. There exists much 
evidence that many real networks possess a scale-free (SF) degree distribution characterized 
by a power law tail given by P(k) ~ k~ x , where k max > k > k min is the degree of a node, 
k m ax is the maximum degree, k min is the minimum degree, and A measures the broadness of 
the distribution 8]. Almost all the studies on networks regarded the links or nodes as iden- 
tical. However, in real networks the links or nodes are not identical but have some "weight." 
As examples the links between computers in the internet network have different capacities 
or bandwidths, resistor networks can have different values of resistance (4], and the airline 
network links connecting pairs of cities in direct flights have different numbers of passengers. 
Many theoretical studies have been carried out on weighted networks 4j, 19J. Recently, several 
studies on real networks with weights on the links, such as the world-wide airport networks 



and the Escherichia coli metabolic networks [10|], have shown that the weights are correlated 
with the network topology and this dramatically changes the transport through them [7, 11]. 
For stance, in synchronization problems, which are very important in brain networks Q, 
networks of coupled populations in the synchronization of epidemic outbreaks [13|, and the 
dynamics and fluctuations of task completion landscapes in causally constrained queuing 
networks 14], the weights could have dramatic consequences for the synchronization 
Synchronization problems deal with optimization of the fluctuations of some scalar field h. 
The system will be optimally synchronized when the fluctuations are minimized. The gen- 
eral treatment to analyze the fluctuations of these processes is to map them into a problem 
of non-equilibrium surface growth via an Edwards- Wilkinson (EW) process on the corre- 
sponding network [l^]. Given a scalar field h on the nodes, that represents the interface 
height at each node, the fluctuations are characterized by the average roughness W(t) of 
the interface at time t, given by W = W{t) = { 1/N Y,?= x % - (h)) 2 } 1/2 , where h { = h^t) 
is the height of node i at time t, (h) is the mean value on the network, N is the system size, 
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and {.} denotes an average over configurations. The EW process on networks is given by 

dhi N 
i=i 

where C^- = Wy is a symmetric coupling strength, {Aij} is the adjacency matrix (Aij = 1 
if z and j are connected and zero otherwise), Wij is the weight on the edge connecting % 
and j, and T]i(t) is a Gaussian uncorrelated noise with zero mean and covariance {T]iT]j} = 
2D5ij5(t — t). Here D is the diffusion coefficient and is taken in general as a constant. For 
non-weighted networks Wy = v = const and thus Eq. ([1]) reduces to the unweighted EW 
equation on a graph given by dhi/dt = v YSj=\ Aij{hj — hi) +7^. Inspired by the results found 
for real networks where the weights are correlated with the topology, Korniss [7| studied 
synchronization for EW processes [see Eq. (CQ)] on SF networks where Wij = {kikjf and fcj 
and kj are the degrees of the nodes connected by a link. Using a mean-field approximation, 
he found that, subject to a fixed total edge cost, synchronization is optimal when (3 = —1, 
and at that point the performance is equivalent to that of the complete graph with the same 
edge cost. Pastore y Piontti et. al [l| used a discrete growth model with surface relaxation 
to the minimum (SRM) in SF networks, which mimics the fluctuation in the task-completion 
landscapes in certain distributed parallel schemes on computer networks, because it balances 
the load. They found that in SF networks with A < 3 the saturation regime of W = W s 
has a logarithmic divergence with N that cannot be explained with the unweighted EW 
equation in graphs, even though in Euclidean lattices the SRM model belongs to the same 
universality class as the EW equation [161 ] . 

In order to understand this discrepancy, in this paper we derive analytically the evolu- 
tion equation for the SRM in random unweighted networks [1] and find that the dynamics 
introduces "weights" on the links. With our evolution equation, which contains non-linear 
terms in the height differences, we recover the logarithmic divergency of W s with N found 
in [l| for SF networks with A < 3. Let us first briefly recall the SRM discrete model 16], 



n 

studied for SF networks by Pastore y Piontti et. al [1]. In this model, at each time step 
a node i is chosen with probability 1/N. If we denote by v% the nearest-neighbor nodes of 
% and j G Vi, then (1) if hi < hj Wj G Vi hi = hi + 1, else (2) if hj < h n Vn ^ j 6 Vi 
=>■ hj = hj + 1. Next we derive the analytical evolution equation for the local height of 
the SRM model in random graphs. The procedure chosen here is based on a coarse-grained 
(CG) version of the discrete Langevin equations obtained from a Kramers- Moyal expansion 



of the master equation 
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19] . The discrete Langevin equation for the evolution of the 



height in any growth model is given by 
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19] 



dht 1 



Gi + rji 



(2) 



at t 

where Gi represents the deterministic growth rules that cause evolution of the node i, r = 
NSt is the mean time to grow a la yer of the interface, and rji is a Gaussian noise with zero 



mean and covariance given by 
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19] 



Mt)v j (t')} = -G l 5 lj 8(t-t / ) 



(3) 



We can write Gi more explicitly as 



A' 



Gi = uJi + Y. A 



(4) 



where u>i is the growth contribution by deposition on node % and Uj is the growth contribution 
to node % by relaxation from any of its j neighbors with 



3&i 



[i _ e(hi - hj)} n [i - ®(hi - h n )\ . 



Here, B is the Heaviside function given by Q(x) = 1 if x > and zero otherwise, with 
x = ht — h s = Ah. Without lost of generality, we take r = 1 and assume that the initial 
configuration of {hi} is random. 

In the CG version Ah — > 0; thus after expanding an analytical representation of Q(x) in 
Taylor series around x = to second order in x, we obtain 
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where c , Ci, and c 2 are the first three coefficients of the expansion of Q(x), Q(kj] 
is the weight on the link ij introduced by the dynamic process, and 
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with C ij — Aij^l^kj^j cLiid Tijfi — j\-ij j^-jyi^l^kj^j . 

In our equation the non-linear terms in the difference of heights arise as a consequence 
of the lack of a geometrical direction and the heterogeneity of the underlying network. This 
result is very different from the one found in Euclidean lattices, where for the SRM model 
the non-linear terms disappear due to the symmetry of the process and the homogeneity of 
the lattice. 

For the noise correlation [see Eq. Q], up to zero order in Ah 

HQ 

we obtain 

{Vi(t)Vj(t')} = 2D{k i )8 ij 8(t - t') with 



1, * 



(7) 



Notice that all the coefficients of the equation depend on the connectivity of node i, i.e., 
on the network topology of the underlying network. This dependence on the topology can 
be thought of as a weight on the links of the unweighted underlying network that appears 
only due to the dynamics on the heterogeneous network. 

Interestingly, the linear terms are different from the EW process as shown below. Keeping 
only the linear terms in Eq. (jSJ), we numerically integrate our evolution equation in a SF 
network using the Euler method with the representation of the Heaviside function given 
by Q( x ) = {1 + tanh[[/(x + z)]}/2, where U is the width and z = 1/2 19]. With this 
representation c = [1 + tanh(C//2)]/2, c x = [1 - tanh 2 ([//2)] U/2, and c 2 = [- tanh(77/2) + 
tanh 3 ([//2)] U 2 /2. We build the network using the Molloy-Reed (MR) algorithm [2o|. In 
Fig. [U we plot W 2 as a function of t, obtained from the integration of Eq. (J2J) using only the 
linear terms of Eq. (jSJ) with D{kj) given by Eq. ((7|) for A = 3.5 and 2.5 and different values 
of iV with k min = 2 in order to ensure that the network is fully connected. For the time step 



integration we chose At <C l/k max according to Ref. [21|. In contrast to the results obtained 
for the EW process [If, W s increases with A" until it reaches a constant value. As shown 
below, this dependence of W s on A" is due to finite-size effects due to the MR construction. 

Now we apply a mean-field approximation to the linear terms of Eq. (jSJ). In this 
approximation we consider 1 <C k min <C k max and disregard the fluctuations. Then 
Ef=iAA/ki « Ef=iaA/a « and Ef=iEn=i,n^nh n /T t « (/*)■ Multiplying 
and dividing Eq. (jSJ) by we can approximate by C^fcj) « fcj J^™^ P(k\ki) Cl(k) dk |7|, 
where P(/c|/cj) is the probability that a node with degree h is connected to another with 
degree k. For uncorrelated networks, P(k\ki) = kP(k)/(k) |8|] does not depend on ki\ then 
Ci(ki) « Ji fcj/(A;) with ii = P(fc) fc fi(fc) dA; . Making the same assumption for T u 
we obtain T^fcj) « 7 2 with J 2 = P(fc) fc (k - 1) fi(A;) rffc. Then the linearized 

evolution equation for the heights can be written as 

dh- 

= Fi{h) + Vi(h) ((h) - hi) + Vl , (8) 

where Fi(ki) = c^+ki l\j (k) represents a local driving force, Vi(kj) = (c\Cq 1 ~ 1 + b)ki is a local 
superficial tension-like coefficient with b = C\(I\ + l<i)j (k), and rji is a Gaussian noise with 
covariance D(ki) = Pj(fcj)/2. This approximation shows the full topology of the network 
through P(k). 

Taking the average over the network in Eq. ©, d(h)/dt = l/NJ2Zi F i = then 
(h) = Ft is linear with t. The solution of Eq. (jSJ) [17] is given by 

ht(t) = f e~ Vi ^ (Fi + Vi(h(s)) +Vi(s)) ds 
Jo 

- (^r) - (^r) ^ + <"> + /„' w * • < 9 > 

Using Eq. (jUJ), the two-point correlation function for t > max {1/Vj} ~ l/k min , is 

{(Hti)- (hmhjfo)- (/»»} = ( F i- F \( F > 

rt 2 rti 



+ / 2 / 1 e -"i(*i--i) e -"i(te-*) {^(s 1 )77 i (s 2 )} d Sl ds 2 . 
Then W s can be written as 

^^Ef^r^E^. do) 
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For SF networks it can be shown that I\, I2 ~ const. + k max exp(—k max const.), where k max ~ 
jyi/C-^- 1 ) f or MR, networks; thus finite-size effects due to the cutoff on these quantities can 
be disregarded. Replacing in the last equation D{ki) by Fj(&j)/2, we obtain 



Wt 



l-2{k) (\) + {k)* (1> 



, v ., 2/ + const. . (11) 

0. 



Notice that, if Fi — 0, D = const, and ^ oc fcj, we recover the EW equation found in 
Using the corrections due to finite-size effects introduced by k max [l| in Eq. (TTTT) . 
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where PV s 2 (oo) = W s 2 (iV — > 00) and q 1 and g 2 are constants. In the inset of Figs. [T] (a) and 
(6) we plot Wg as function on N and the fitting obtained from Eq. ( f!2l) . The agreement with 
the scaling form, Eq. ( fl2|) . is excellent. Thus, the linear approximation can only explain the 
finite-size effects due to the MR construction but fails to predict the logarithmic divergency 
of W a with iV for A < 3 found in Ref. [l|. Next we show that the non-linear terms are 
responsible for this behavior. We integrate our evolution equation for SF networks with the 
linear terms and only the first non-linear term [see Eq. ([5])] due to the numerical instability 
produced when we try to incorporate all of them. Even with only one non-linear term, we 
recover the logarithmic divergency of W s with N for A < 3. The results of the integration 
are shown in Fig. [2j where we plot W as a function of t for (a) A = 3.5 and (b) A = 2.5 and 
different values of N. In the inset figures we plot W s as a function of N. We can see that, 
for A = 3.5, W s increases but asymptotically goes to a constant and all the N dependence 
is due to finite-size effects. However, for A = 2.5 we found a logarithmic divergency of W s 
with N [1], as shown in the inset of Fig. [2] (b), where we plot W s as a function of iV on a 
log-linear scale. The fit of W s with a logarithmic function for A = 2.5 shows the agreement 
between our results and those obtained for the SRM model in SF networks for A < 3. 
Discrepancies between behaviors in regular Euclidean lattices and Euclidean lattices after 
addition of random links were found before in {22]. 

In summary, we derived analytically the evolution equation for the SRM model and found, 
surprisingly, that even when the underlying network is unweighted the dynamics introduces 
weights on the links that depend on the topology. We also found that the linear terms can 
explain only finite-size effects due to the MR construction. The linear mean-field approx- 
imation shows clearly the effects of the topology on the dynamics and the corrections due 



to finite-size effects. When non-linear terms on SF networks are considered, new numerical 
integration algorithms are needed in order to avoid numerical instabilities. This is still an 
open problem to be solved in the future. With all the linear terms and one non-linear term, 
we recovered the logarithmic divergency of W s with iV of the SRM model for A < 3. Our 
analytic procedure can be also applied to any other growth model. 

We thank A. L. Pastore y Piontti for useful discussions and comments. This work has 
been supported by UNMdP and FONCyT (Pict 2005/32353). 
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FIG. 1: W 2 as a function of t from the integration of the evolution equation using the linear terms 
for N = 256 (o), 384 (□), 512 (o), 768 (A), and 1024 (y). A = (a) 3.5 and (6) 2.5. In the inset 
figure we plot W 2 vs N in symbols. The dashed lines represent the fitting with Eq. (|12|) . obtained 
by considering the finite-size effects introduced by the MR construction. For all the integrations 
we used U = 0.5 and typically 10 000 realizations of networks. 
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FIG. 2: Wasa function of t from the integration of the evolution equation using the linear terms 
and the first non-linear term for N = 384 (o), 512 (□), 768 (o), 1024 (A), and 1536 (y). A = (a) 3.5 
and (b) 2.5. In the inset figure we plot W s vs N in symbols. The dashed lines represent the fitting 
with Eq. CE]) in (a) and W s ~ hi TV in (b). 
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